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ABSTRACT 

We examine masses of hosting haloes of two photometricahy-selected high- 2; galaxy 
samples: the old passively-evolving galaxies (OPEGs) at z 1 and Lyman Break 
Galaxies (LBGs) at z ~ 4 both taken from the Subaru/XMM-Newton Deep Survey 
(SXDS). The large survey area of the SXDS (1 deg^) allows us to measure the angu- 
lar two-point correlation functions to a wide separation of > 10 arcmin with a good 
statistical quality. We utilize the halo model prescription for estimating character- 
istic masses of hosting haloes from the measured large-scale clustering amplitudes. 
It is found that the hosting halo mass positively correlates with the luminosity of 
galaxies. Then, adopting the extended Press-Schechter model (EPS), we compute the 
predictions for the mass evolution of the hosting haloes in the framework of the cold 
dark matter (CDM) cosmology in order to make an evolutionary link between the 
two galaxy samples at different redshifts and to identify their present-day descendants 
by letting their haloes evolve forward in time. It is found that, in the view of the 
mass evolution of hosting haloes in the CDM model, bright (i' < i\ -I- 1) LBGs are 
consistent with being the progenitor of the OPEGs, whereas it is less likely that the 
LBG population, as a whole, have evolved into the OPEC population. It is also found 
that the present-day descendants of both the bright LBGs and OPEGs are likely to be 
located in massive systems such as groups of galaxies or clusters of galaxies. Finally, 
we estimate the hosting halo mass of local early- type galaxy samples from the 2dF and 
SDSS based on the halo model and it turns out that their expected characteristic mass 
of hosting haloes is in good agreement with the EPS predictions for the descendant's 
mass of both the bright LBGs and OPEGs. 

Key words: galaxies: formation — galaxies: haloes — galaxies: high-redshift — 
cosmology: theory — dark matter 



1 INTRODUCTION 

Galaxies consist of two major ingredients, namely, the bary- 
onic and dark matter. There are pieces of observational evi- 
dence that properties and evolution of the baryonic compo- 
nent are not independent of those of the dark matter, but 
are closely correlated. For instance, early/late type galax- 
ies dominate the bright/faint part of the galaxy luminos- 
ity function, which suggests that the star formation his- 



* Based on data collected at Subaru Telescope, which is operated 
by the National Astronomical Observatory of Japan. 



tory and/or morphology are related to the total mass of 
the system, provided that the luminosity is basically propor- 
tional to the total mass. Another example is the morphology- 
density relation (Dressier 1980; Postman & Geller 1984), 
which may indicate the possible influence of the mass of 
hosting halo on the galaxy formation on group/cluster scale 
as well. It follows from these pieces of evidence that the mass 
of the hosting halo could be one of the most fundamental 
quantities in the galaxy formation. Therefore, in order to 
understand the galaxy formation, the evolution of baryonic 
and dark matter should not be treated as independent pro- 
cesses but should be investigated in a unified way. 
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The baryonic contents (stars and gas) of galaxies have 
been weU studied by both photometric and spectroscopic ob- 
servations and have been extensively analyzed with a help 
of well developed theoretical tools such as stellar popula- 
tion synthesis models (e.g., Bruzual & Chariot 2003). On 
the other hand, the dark matter component of hosting halos 
has been less studied observationally due to the technical 
difflculties in probing it since it requires dynamical tracers 
or gravitational lensing. Theoretical models for the evolu- 
tion of dark matter, however, have been well developed with 
numerical simulations and analytical models such as Press- 
Schechter prescriptions (Press & Schechter 1974). 

It has been known that a large-scale clustering ampli- 
tude of galaxies provides a unique way of estimating their 
hosting halo mass in a statistical manner. This makes use of 
a finding from theoretical studies of dark-matter structure 
formation in the cold dark matter (CDM) model that a clus- 
tering amplitude of dark matter haloes depends monotoni- 
cally and strongly on the halo mass (Mo & White 1996). By 
measuring the clustering amplitude of galaxy populations 
selected from a wide-field survey, we can relate properties of 
the populations such as stellar mass and/or star formation 
rate with the mass of their hosting halo (e.g., Giavalisco & 
Dickinson 2001; Adelberger et al. 2005). It, in turn, enables 
us to explore an ancestor-descendant connection of different 
galaxy populations at various redshifts from the viewpoints 
of both the star formation and the dark matter assembly 
(Moustakas & Somerville 2002; Hamana et al. 2004; Ouchi 
et al 2004b). 

This approach is exactly the one we take, in what fol- 
lows, to explore the evolution of the old passively-evolving 
galaxies (OPEGs) at z ~ 1 and Lyman Break Galaxies 
(LBGs) at z ~ 4. The OPEGs are selected by optical colour 
criteria (see >13.1II being consistent with a passively evolv- 
ing galaxy with the major star-formation epoch at z > 2. 
Since the age of the universe at z = 2 is only ~ 2.3/i~^Gyrs, 
progenitors of the OPEGs must have experienced an active 
star formation phase at some high redshifts of 2 > 2. The 
LBGs, which are in general actively star forming galaxies at 
high redshifts, are naturally a strong candidate of a major 
population of the progenitor of such OPEGs. While there is 
a fraction of dusty star-forming galaxies or red old galax- 
ies with little star formation, which are not selected by the 
LBG criteria at z ~2-2.5 (e.g., Franx et al. 2003; Reddy 
et al. 2005), LBGs are still most populous among them. In 
addition, the OPEGs are considered as strong candidates 
for progenitors of the present-day early-type galaxies since 
the simple passive evolution in their stellar populations can 
easily link between them. We examine relations between the 
two galaxy populations in terms of the evolution of hosting 
dark matter haloes. Our primary question addressed in this 
paper is "supposing the LBGs are a major progenitor of the 
OPEGs, are their hosting halo masses compatible with the 
mass evolution of dark matter haloes in the CDM structure 
formation model ?" Also, we discuss possible present-day 
descendants of those high-z galaxies from the viewpoint of 
the masses of hosting haloes. 

We determine the clustering amplitude at large scales 
with angular two-point correlation functions of the OPEGs 
and LBGs measured, with a good accuracy, from a wide field 
(1 deg^) deep imaging data-set of Subaru/XMM-Newton 
Deep Survey (SXDS). Multi-band colour selection tech- 



niques successfully isolate those galaxies with well calibrated 
redshift selection functions, which enable us to compute an 
accurate prediction for the corresponding dark matter an- 
gular correlation function via the Limber's projection. This 
large data-set and selection techniques allow us to compute 
a large-scale galaxy bias in a robust manner. Comparing 
the measured large-scale bias with the ones predicted by 
the halo model, we place a limit on the hosting halo mass 
of galaxy populations. Here we adopt an empirically pa- 
rameterized model for the halo occupation function (HOF) 
which describes statistical relations between galaxies and 
their hosting haloes. Then, utilizing the extended Press- 
Schechter (EPS) prescriptions, we examine the evolution of 
halo mass in the framework of the CDM model. In this way, 
we compare hosting halo masses of each galaxy population 
at different redshift, and explore their connection. 

The outline of this paper is as follows. Section 2 de- 
scribes models and basic equations. Section 3 summarizes 
observational data that are used to place a constraint on 
the hosting halo mass. Results are presented in section 4. 
Finally, section 5 is devoted to a summary and discussion. 

Throughout this paper, we adopt a flat ACDM cos- 
mology with the matter density ilm = 0.3, the cosmolog- 
ical constant JIa ~ 0.7, the Hubble constant Ho = lOO/i 
km s~^ Mpc~^ with h — 0.7, and the normalization of the 
matter power spectrum erg = 0.9. We adopt the fitting func- 
tion of the CDM power spectrum of Bardeen et al. (1986). 
We present magnitudes in the AB system. 



2 MODELS 

2.1 Dark matter angular correlation function 

We quantify a clustering amplitude of a population of galax- 
ies by comparing their angular two-point correlation func- 
tion with the corresponding dark matter correlation func- 
tion. Let q{z) be a normalized redshift selection function of 
a population of galaxies being considered, the dark matter 
angular two-point correlation function is computed from the 
dark matter power spectrum (Pdm) via the Limber projec- 
tion (see e.g, chapter 2 of Bartelmann & Schneider 2001): 

LODMie) = J drq\r) J PuM{k,r) Jo[fK{r)9k], (1) 

where Jo (a;) is the zeroth-order Bessel function of the first 
kind. For the spatially flat cosmology (fim -I- SIa ~ 1) as we 
consider throughout the present paper, the radial function 
fKir) is equivalent to r, and r = r{z) is the radial comoving 
distance given by r{z) = c/Ho dz' [^^{1 + z'f + nA]~^^'^ ■ 
We use the nonlinear fitting function of the CDM power 
spectrum by Peacock & Dodds (1996). 

2.2 Halo model 

We utilize the halo model approach for estimating a charac- 
teristic mass of hosting haloes from the measured clustering 
amplitude of galaxies. Here, we summarize several expres- 
sions which are most relevant to the current analysis. See 
Hamana et al. (2004; and references therein) for details of 
the halo model. 

We adopt a simple parametric form for the average 
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number of a given galaxy population as a function of the 
hosting halo mass: 



(M > M™„) 
(M < M™„) 



(2) 



which is characterized by the three parameters: The min- 
imum mass of haloes which host the population of galax- 
ies (Mmm), a normalization parameter which can be inter- 
preted as the critical mass above which haloes typically host 
more than one galaxy (Mi), and the power- law index of the 
mass dependence of the galaxy occupation number (a). For 
the dependences of the HOP parameters on the shape of the 
two-point correlation function, see Hamana et al. (2004; and 
references therein). 

We introduce the average mass of hosting halo 
(weighted by the number of member galaxy) by 



(A/host) = 



JdM M iVg(M)nhaio(M) 



(3) 



where nhaio(Af) denotes the halo mass function for which 
we adopt the fitting function of Sheth & Tormen (1999) 

Since the halo model assumes the linear halo bias (Mo & 
White 1996) and Ng{M) solely depends on the halo mass, 
the galaxy bias on large scales (the scales larger than the 
virial radius of hosting haloes) is given by the galaxy number 
weighted halo bias: 



Og,L 



J dM &haio(M)jVg(M)nhaio(M) 
/ dM iVg(M)nhaio(M) 



(4) 



where bhaio(-M) is the halo bias, for which we adopt the 
fitting function of Sheth & Tormen (1999). Note that bg^L 
does not depend on Mi but only on Mmin and a. 

Since galaxies considered in the following sections are 
distributed over a redshift interval, we take into account the 
redshift evolution of a quantity X{z) (represents for or 
(Mhost)) by computing its redshift average: 

/ dz [dV/dz] q{zfX{z) 



X 



J dz [dV/dz] q{z) 



(5) 



where k = 1 for (Mhost) while fe = 2 for 6g,L! and dV/dz 
denotes the comoving volume element per unit solid angle: 
dV/dz = c/Hq [rim(l + z^ + VL\\^^/^, again for the spa- 
tially flat cosmology. 



2.3 Extended Press-Schechter model 

The extended Press-Schechter (EPS) formalism was devel- 
oped by Bond et al. (1991), Bower (1991) and Lacey & Cole 
(1993). Since the EPS model provides a way to treat halo 
mergers, which play an important role in the structure for- 
mation in the hierarchical CDM model, it has been widely 
applied to analytical and semi-analytical studies of the struc- 
ture formation. We utilize the EPS model for making a sta- 
tistical estimate of growth of halo masses. A key expression 
for this is the conditional probability P2[Mt2, Z2\Mti, z\) 
that a material in a halo of mass Mti at zi will be in 
a halo of mass Mt2 {Mt2 > Mti) at a later redshift 22, 
leading to an expression for the conditional mass function 
n2{Mt2,Z2\Mti,zi). 

Here we summarize only expressions which are directly 



relevant to our analysis, see above references for their deriva- 
tion. Let Sc and a be the critical density threshold for a 
spherical perturbation to collapse and the RMS density fluc- 
tuation smoothed over a region enclosing a mass M, respec- 
tively, the conditional probability is 



P2{Mt2,Z2\Mtl,Zl)dMt2 
_ 1 Sc2 {Scl — Sc2) 

(cr|(5cl 



27r 



cr2(o 



X exp 



■ (7idc2 ) 



dMt2 



3/2 



dMt: 



(6) 



where the subscripts 1 and 2 stand for epochs z\ and 22, 
respectively. The conditional mass function is given from 
this by 

n2(Mt2,22|Mti,2l)dMt2 OC rr^P2(Mt2,22|Mtl,Zl)dMt2.(7) 

Mt2 



3 DATA 

We use two samples of photometrically-selected galaxies 
from B, R, i' and z' imaging data of the Subaru/XMM- 
Newton Deep Survey^ (SXDS). The limiting magnitudes for 
the 3(7 detection of an object in a 2 arcsec diameter aperture 
are B ~ 28.3, R ~ 27.6, i' ~ 27.5 and 2' ~ 26.5 (Furusawa 
et al. in preparation). The seeing size of those images is 
~ 0.8". We avoid the edges of the field and the area affected 
by bright sources. The final effective area used in this paper 
is ~ 1 deg^. 

Below we summarize basic properties of the two galaxy 
samples, and present angular two-point correlation functions 
of those galaxies. 



3.1 Old passively-evolving galaxies at 2 ~ 1 

The photometric properties of the old passively-evolving 
galaxies (OPEGs) in the SXDS fields are presented in Ko- 
dama et al. (2004) and Yamada et al. (2005). For the sample 
selection of OPEGs, we follow the definition of Yamada et 
al. (2005) which imposes the two criteria: 0.8 < i' — 2' < 1.2 
and -0.052' -1-3.01 < R- z' < -0.032'-|-2.49, on the 2'-band 
selected catalog. This criteria effectively isolate the galaxies 
with star- formation epoch greater than 2/ > 2 and located 
at 0.9 < z < 1.1. This photometrically-selected sample is 
composed of 4,118 OPEG candidates with 2' < 25.0 dis- 
tributed in the effective area of 1.03 deg^. It is found from 
the spectroscopic observations for the 93 OPEG candidates 
with 19 < 2' < 22 that the 73 objects he between 2 = 0.87 
and 1.12, and the 4 objects lie outside of the redshift interval, 
and the remaining 16 objects have undetermined redshift be- 
cause of no usable feature in the spectrum or a poor S/N. 
Figure 8 of Yamada et al. (2005) shows the redshift distribu- 
tion function of OPEGs obtained from the spectroscopically- 
identified objects. The contamination fraction is estimated 
to be between 0.05 (=4/[73-f4]) and 0.22 (=20/93). We as- 
sume the contamination fraction of fc =0.1 in the follow- 
ing analyses. The luminosity function is well fitted by the 



^ See http:/ / www. naoj .org/Science / SubaruPro ject /SDS/ [ for de- 
tails of the SXDS project. 
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Figure 1. Upper panel: Filled circles with error bars show the an- 
gular two-point correlation function for OPEGs with z' < 23.8(= 
+2). Note that the plotted correlation function has been cor- 
rected for the contamination (see text). The dotted line shows 
the CDM model prediction for the dark matter angular two-point 
correlation function computed via equation 0, where the non- 
linear fitting function for the CDM power spectrum by Peacock 
& Dodds (1996) is used. Lower panel: The corresponding galaxy 
bias defined by equation jgj is plotted. The horizontal solid line 
shows the large-scale bias factor computed averaging the bias over 
an interval l' < 9 < 10', and the dashed lines show its l-cr error. 



Schechter function with = (4.26 ±0.42) x 10~^/i^Mpc"^, 
aLF = -0.67 ± 0.07 and Ms, = -21.38 ± 0.10 (Yamada et 
al. 2005). The latter magnitude corresponds to zi — 21.8 at 
z = l. 

The angular two-point correlation functions are com- 
puted using the pair-count estimator formulated by Landy 
& Szalay (1993): iUg{e) = [DD{e)-2DR{e)+RR{e)]/RR{e). 
In so doing, we distribute the same number of random sam- 
ples with the same geometrical constraint as of the data 
sample. We repeat 100 random re-samplings, and the mean 
and RMS among the 100 measurements are taken as the 
mean correlation signal and l-cr error. The effect of the con- 
tamination on the two-point correlation function is corrected 
by multiplying by a factor of 1/(1 — /c)^, where the contam- 
inants are assumed to be randomly distributed. A measured 
angular correlation function is plotted in the top panel of 
Figure together with the dark- matter angular correlation 
function computed with the same redshift selection function. 
Note that the integral constraint (Groth & Peebles 1977) is 
estimated to be ~ 7 x 10~^ for a single power-law model 
[L0g{6) oc 6^] with the power-law index of /3 = — 1, we may 
thus ignore an effect of a finite field size as far as 9 < 10' is 
concerned. 

We define the galaxy bias by 




(8) 



23 24 
z' [AB mag] 

Figure 2. Lower panel: The large-scale bias factor (see test for 
the definition) for the OPEGs as a function of limiting magni- 
tude. Upper panel: (Mhost) computed from the halo model (with 
a fixed a = 0.8) for the corresponding large-scale bias factor (see 
j4]for details). The upper abscissa axis indicates the magnitude 
difference from the z' value. 



and is plotted in the bottom panel of Figure As is 
shown there, on scales below ~ 1 arcmin, the bias decreases 
with the separation, while on larger scales it flattens. The 
comoving angular length of the transition scale (~l') is 
~ 0.7/i^^Mpc (at z = 1) which corresponds to the virial 
radius of the halo with the mass ~ 2 x IQ^'^ Mq at 2 = 1. 
This mass coincides with the characteristic mass of the host- 
ing halo ((Mhost)) predicted by the halo model analysis in 
the next section. It may follow from this that the shape of 
the bias function is basically understood by the standard 
halo model picture: The small-scale clustering arises from 
galaxy pairs located in the same halo, while the large-scale 
clustering arises from galaxy pairs located in two different 
haloes. 

Since the measured bias flattens on large-scales as ex- 
pected by the halo model, we define the large-scale bias fac- 
tor (6l) as an averaged bias over 1' < 9 < 10' (illustrated in 
the bottom panel of Figure Ql. We compute the large-scale 
bias factor for OPEG samples with different limiting mag- 
nitudes and show in the bottom panel of Figure H Broadly 
speaking, the large-scale bias factor is a decreasing function 
of the magnitude, though the significance is not high due to 
large error bars. 

Let us devote, in passing, a little space to compare our 
measurement with previous studies. We derive the comov- 
ing correlation length, ro, from the power-law fitting model 
of the angular two-point correlation function with the Lim- 
ber equation (Peebles 1980). The correlation length ro is the 
normalization of the spatial two-point correlation function, 
^ = (r/ro)~^, where 7 is related by 7 = /? -f- 1. We find 
ro = 4.7 ± 0.3/i"^Mpc and 5.7 ± 0.2/i"^Mpc with the best- 
fit beta {j3 = 1.1) and the fixed beta (/3 = 0.8), respectively. 
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Figure 3. Same as Figure^Jbut for the LBGs with i' < 27.0(= 
<+2). 




26 
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Figure 4. Same as Figure 1^ but for the LBGs. 



This is consistent with ro — 5.0±0.3ft^ Mpc obtained from 
{R — I) colour selected red galaxies (7 — 18 — 24) at z ~ 0.85 
most of which have early-type spectra (Coil et al. 2004). In 
contrast, it is much smaller than ro = (8 — 13)^~^Mpc ob- 
tained from (optical— NIR) colour selected "Extremely Red 
Objects" (EROs) at 1 < 2 < 2 (Daddi et al 2001; Mc- 
Carthy et al. 2001; Firth et al. 2002; Roche et al. 2002; 
Miyazaki et al. 2003) . The luminosity segregation of cluster- 
ing may not explain the difference of ro, since, if we assume 
an early-type spectrum, the limiting magnitudes of their 
ERO searches reach to a comparable depth to our OPEG 
sample {z' < 23.8). The reason for the significant difference 
in ro is not clear. Coil et al. (2004) suggest that their red 
galaxies selected in optical bands may be a less-extreme ver- 
sion of EROs. Another possible reason would be an intrinsic 
evolution in the clustering strength due to the difference in 
redshift ranges of the samples. While our sample and Coil 
et al.'s (2004) sample are both limited to the small redshift 
range at 2 ~ 1, other "R — K" or "J — K" EROs sam- 
ples have much wider range in redshift up to 2: ~ 2 The 
clustering strength of the red galaxies can evolve strongly 
between z ~ 2 and 1. Furthermore, the SXDS survey probes 
much larger volume than the other surveys and the effect of 
field-to-field variation is expected to relatively small. Fur- 
ther investigating into the origin of the difference is, how- 
ever, beyond the scope of this paper, and we leave it for a 
future work. 



3.2 Lyman break galaxies at ^ ^ 4 

The sample selection and clustering properties of the Lyman 
break galaxies (LBGs) at z ~ 4 are described in Ouchi et al. 
(2005). Briefiy, the LBGs are selected from i'-band selected 
catalog by three criteria: B — R > 1.2, R — i' < 0.7 and 
B — R > 1.6{R — i') + 1.9. The number counts of this sam- 



ple is given in Table 1 of Ouchi et al. (2005). The redshift 
distribution function has the mean of 2; ~ 4 and width of 
Az ~ 0.5 as is shown in the top panels of Figure 12 of Ouchi 
et al. (2004a). It is found from the spectroscopic follow-up 
observations for the 63 photometrically-selected LBGs that 
the 60 objects lie between z — 3.5 and 4.5 (Ouchi et al. 
2005 and references therein). The contamination fraction is 
thus estimated to be fc = 0.05(= [63 — 60]/63). Ouchi et 
al. (2004a) computed the luminosity function of the LBGs 
selected from the Subaru Deep field with the same colour 
selection criteria and have found i» = 25.0 ± 0.1. 

The angular two-point correlation functions of the 
LBGs are computed in the same procedure as was done for 
the OPEGs and are plotted in the top panel of Figure|3 The 
integral constraint is estimated to be ~ 7 x 10~^ for a sin- 
gle power-law model with the power-law index of /3 = —0.8, 
we may thus ignore an effect of a finite field size as far as 
9 < 5' is concerned. The galaxy bias defined by equation JSJ 
is plotted in the bottom panel. A transition from the small- 
scale decreasing part to the large-scale flat part is found at 
the scale of ~10 arcsec (see Ouchi et al. 2005 for further de- 
tail discussions on the shape of the bias function. Note that 
they do not correct for the effect of the contamination be- 
cause of unknown clustering property of the contaminants, 
for which we assume the random distribution. As a con- 
sequence, their measured biases are smaller than ours by 
about 5 percents). This is transformed to the comoving an- 
gular scale of ~0.2/i~^Mpc at 2; = 4, which corresponds to 
the virial radius of the halo with the mass ~ 1 x 
This mass, again, coincides with the characteristic mass of 
hosting haloes ((Mhost)) predicted by the halo model anal- 
ysis in the next section. Thus the behavior of the bias is 
understood by the standard halo model picture. 

The large-scale bias factor is computed in the same 
manner as was done for the OPEGs except for using a dif- 
ferent separation range of 0.5' < d < 5' (illustrated in the 
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Figure 5. Halo model predictions for OPEGs. Lower panel: The 
large-scale bias bg i, defined by equation jlj for three values of a 
(1.0, 0.8 and from upper to lower). The horizontal lines show 
the measured large-scale bias factor for the OPEG sample with 
z' < 23.8(= z'^ + 2), the solid and dashed line represent the 
mean and l-tr range, respectively. Upper panel: The mean mass 
of hosting haloes defined by equation for a = 1.0, 0.8 and 
from upper to lower. 



bottom panel of Figure 1^. This bias factor is calculated as a 
function of limiting magnitudes and is plotted in the bottom 
panel of Figure 0] As was reported by Ouchi et al. (2005), 
a trend that the large-scale bias factor decreases with the 
luminosity is observed. 



4 RESULTS 

We define fiducial subsamples of the two galaxy populations 
for the later analyses by setting the magnitude limit to 2 
magnitude fainter than the m« value [z'^, — 21.8 for the 
OPEGs and = 25.0 for the LBGs). This is chosen so that 
it is well deeper than the m« value and at the same time 
it is well brighter than the completeness limit of our imag- 
ing data {z' ~ 26.5 for the OPEGs and i' ~ 27.5 for the 
LBGs). The former is imposed so that the samples do not 
contain only the brightest objects, which would not be rep- 
resentative of the whole population. The latter is imposed to 
reduce contaminations due to errors in the luminosities and 
colours. We note that the number densities of the subsam- 
ples are nopEG (2' < 23.8) = (4.7±0.47)xl0"3/i3Mpc-3 and 
nLBG(i' < 27.0) = (1.1 ±0.2) x lO'^/i^Mpc-^ for OPEGs 
and LBGs, respectively. Therefore with this selection, the 
LBGs are as 2.3 times numerous as the OPEGs. It should 
be noted that statistical properties (e.g., the number density 
and clustering amplitude) of a galaxy sample would be sen- 
sitive to the selection criteria (e.g., the magnitude limit and 
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Figure 6. Same as Figure 151 but for the LBGs. Measured large- 
scale bias factor and number density shown by the horizontal lines 
are taken from the sample i' < 27.0(= + 2). 



Table 1. Constraints on (A/jjost) obtained from the halo model 
analysis. Unit of the mass is in lh~^MQ]. values of (Afhost) 
estimated from the mean of b^. (2) the interval of (Mhost) from 
the l-cr range of fej,. 
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colour selection) . Therefore, when one attempts to compare 
properties of two (or more) galaxy samples, it is necessary 
to properly define samples of galaxies with a quantity es- 
sential to their nature. In our case, the luminosity would 
be the most relevant quantity among (a few) controllable 
parameters we have in hand. Since the definition of our sub- 
samples is based on a somewhat arbitrary magnitude limit, 
we shall look into an effect of this definition by changing the 
magnitude limit for subsampling. 



4.1 Halo model analyses 

We estimate the hosting halo mass by comparing the mea- 
sured large-scale bias factor with the halo model prediction. 
To do this, we proceed as follows: First, we compute the 
halo model prediction for the large-scale bieis 69, l defined 
by equation @ for a given value of a as a function of Mmin 
(plotted in the lower panel of Figures|3and|SJ . Then, search- 
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ing for an interval of Mmin, where the predicted bg,L and the 
l-cr interval of the measured large-scale bias factor intersect, 
we have a constraint on Mmin- Since for a given a, (Mhost) 
and Mmin have a one-to-one correspondence (see the upper 
panel of Figures]^ and |SJ, the constraint on Mmin is imme- 
diately translated into the constraint on (Mhost)- We take 
three values of a = 1, 0.8 and 0. Here a < 1 is preferred from 
semi-analytic models of the galaxy formation (e.g., Kravtsov 
et al. 2004), as well as from the halo model analysis for LBGs 
(Hamana et al. 2004). We take a = 0.8 as a fiducial value. 
The case q = corresponds to an extreme case where every 
halo with a mass greater than Mmin always has one galaxy. 

The constraints obtained for our fiducial galaxy samples 
[z < 23.8 for the OPEGs and i' < 27.0 for the LBGs) 
are summarized in Table It is important to notice that 
the change in a does not make a significant change in the 
preferred interval of (Mhost) except for the extreme case of 
a = for the OPEGs. Therefore , the constraint is not very 
sensitive to the uncertainty in a. 

Let us look into how the characteristic hosting halo 
mass varies with the limiting magnitude for the sample se- 
lection. The upper panels of Figures |5| and 2] show the con- 
straint on (Mhost) as a function of the limiting magnitude. 
In both galaxy populations, a trend of decreasing (Mhost) 
for a fainter limiting magnitude is observed. These similar- 
looking trends, however, may arise from difi^erent physical 
origins as explained below. First, for the OPEGs, the ob- 
served z'-band corresponds approximately to the rest-frame 
B-band. Although the B-band luminosity is affected by on- 
going/recent star formation, only a weak star formation ac- 
tivity would make the colour of galaxy blue and pushes it 
outside of our OPEG selection criteria (Yamada et al. 2005) . 
Therefore the _B-band luminosity is a good measure of the 
stellar mass of the OPEGs. Accordingly, the observed trend 
is considered as a result of a correlation between the stel- 
lar mass and the hosting halo mass. On the other hand, for 
the LBGs, the observed i'-band corresponds to ~1500A in 
the rest-frame, where the luminosity is most sensitive to the 
star formation. The LBGs are generally in an active star for- 
mation phase, also it was reported that the stellar mass of 
LBGs is poorly correlated with the UV luminosity (Shapley 
et al. 2001). Therefore, the observed luminosity dependence 
of hosting halo mass of LBGs may suggest a presence of 
a correlation between the hosting halo mass and the star 
formation activity. 

4.2 EPS model analyses 

We use the conditional mass function n2{Mt2, Z2\Mti, zi) 
(equation |7J derived from the EPS formalism to predict the 
mass evolution of hosting haloes in the framework of the 
CDM cosmology. For the parameters of the haloes at earlier 
epoch zi, we take the l-cr interval of (Mhost) obtained from 
the halo model analysis with a = 0.8 as a range of Mti, 
and we set zi = 4 and zi = 1 for the LBGs and OPEGs, 
respectively. Then for a certain later epoch Z2, we compute 
?^2(Mt2, Z2\Mti, Zi) as a function of Mt2 (see inserts of Fig- 
ures |7JJ5| for examples) . Assuming that the mass assembly 
history of hosting haloes is not biased toward any specific 
merging path, the conditional mass function can be regarded 
as the probability distribution function (PDF) of the mass 
of descendant haloes A'It2- We define a 68% confidence in- 
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Figure 7. Inserts show arbitrary normalized conditional mass 
functions n2{Mt2, Z2\Mti, zi) for Z2 = 1 (left panel) and (right 
panel), where parameters of the earlier halo is taken from the 
limit obtained by the halo model analysis of the LBGs (with 
a = 0.8): zi = 4, Mti = 1.2 X IO^'^H-^Mq (red), 8.3 X lO^^ 
(blue), and 1.7 X lO'^^ (green). These conditional mass functions 
are considered as the probability distribution function of the mass 
of the descendant. The main plot shows 68% confidence intervals 
of Mt2 as a function of the redshift (22). Three intervals shown by 
coloured hatches are for the same three Mti values as the inserts. 
Filled circles with error bars show the limits of (-Mhost) obtained 
from the halo model analysis (with a = 0.8) for the LBGs and 
OPEGs. 



terval of Mt2 as the EPS model prediction for the mass of 
the descendant halo. 

Such predictions are made for our LBGs and OPEGs 
samples, and are shown in Figures [TJI^I as hatched regions. 
Three tracks correspond to the three different values of Mn 
taken from the mean and upper/lower 1-a values for the 
(A/host) denoted by the symbols with error bars (a — 0.8 
cases in Table 0. The inserts show the PDFs of Mt2 at 
z = (and z — 1). It is observed in the inserts that the 
PDFs have a wide spread, implying a wide variety of the 
mass assembly history of hosting haloes. We note that since 
the PDFs are skewed toward a larger mass, the mean of the 
distribution is greater than the mode. 

Let us look closely at each result. First of all. Figure Q 
compares the expected descendant mass of hosting haloes of 
the LBGs with the (Mhost) of OPEGs. Both galaxy samples 
are selected with our fiducial magnitude limit of m < m*-|-2. 
As is evidently shown, the EPS predictions for the mass of 
the LBG descendant are slightly smaller than the predicted 
mass range of the OPEGs, only the track from the upper 
limit of the LBG haloes is compatible with the OPEGs. 

Is there a selection criterion for LBGs or for OPEGs 
which results in a better agreement? It is found from the 
upper panel of Fig. |21 that choosing a fainter limiting mag- 
nitude for the OPEG sample lowers the (Afhost) very little 
and can hardly solve the incompatibility. On the other hand, 
choosing a brighter limiting magnitude for the LBG sample 
raises the (Mhost) as shown in the upper panel of Figure |3 
Indeed, the galaxy sample selected with i' < 26.0(= i(, -I- 1) 
results in a very compatible halo mass with the (Mhost) of 
the OPEGs (with z' < 23.8) as shown in Figure H There- 
fore, in the view of the mass evolution of hosting haloes, it 
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Figure 8. Same as Figure|7|but the LBGs sample is selected with 
i' < 26.0. The limit of (Mi^ast) obtained from the halo model 
analysis with a = 0.8, which is taken by Mti, is Mti = 2.2 X 
W^^H-^Mq (red), 1.2 X lO^^ (blue), and 3.7 X lO^^ (green). 
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Figure 9. Same as Figure □ but the limit for OPEGs obtained 
from the halo model analysis (with a = 0.8) is taken by param- 
eters of the earlier halo: zi = 1, Mti = 1.2 X lO^^h'^MQ (red), 
6.3 X 10^^ (blue), and 1.9 X 10^* (green). In the insert, Z2 is taken 
by 0. 




may be safely concluded that the bright {i' <«', + !) LBGs 
are consistent with being the progenitor of the OPEGs, 
whereas it seems less likely that the LBGs population, 
as a whole, has evolved into the OPEG population. Note 
that the corresponding number density is computed to be 
nLBG(i' < 26.0) ~ (3.8 ±0.3) x IQ-'^h^Mpc'''^ , closer to the 
number density of nopEG(z' < 23.8). 



4.3 Predictions for the present-day descendants 

Turn next to predictions for the present-day descendants of 
the two galaxy populations plotted in Figures |7JJ^ As is 
shown in the inserts of Figures [7| and |S| the PDF of the 
present-day descendants of the LBGs have a very broad 
spread due to the wide variety of the mass assembly history 
over ~ 8/i~^Gyrs. The predicted mass range for the LBG 
sample with i' < 27.0 is broadly consistent with that of the 
OPEGs, though the PDF of the LBGs extends to a smaller 
mass (say M < 10^^/i"^Mo) where the PDF of the OPEGs 
has little probability. On the other hand, the predicted mass 
range for the LBG sample with i' < 26.0 agrees better with 
that of the OPEG descendants. In this case, both the PDFs 
computed from the central (Mhoat) value (the middle one of 
the three cases in the inserts) peak at ~ 2 x 10^^ Mq, 
which corresponds to the mass scale of groups of galaxies. It 
is important to note that since the PDFs are skewed strongly 
toward a larger mass, a certain fraction of haloes is expected 
to evolved into more massive haloes with Mt2 > IO^/i'^Mq 
(see the inserts of Figures |S] and |UJ. 

It is interesting to compare those EPS predictions for 
the halo masses of the present-day descendants with local 
galaxy samples. Since both the LBGs and OPEGs are fre- 
quently argued as the strong candidates for the progenitor 
of the present-day early-type galaxies, we take two cluster- 
ing analyses of the early-type galaxy samples: One from the 
the 2dF (Madgwick et al. 2003), and the other from the 
SDSS (Zehavi et al. 2002). Madgwick et al. (2003) found 
a^^ ~ 1.1 ± 0.1 for "passive" galaxies, of which the lu- 
minosity range is —16.5 < Mi,j — 51og/i < —22 with 
Mb J,* = -19.7. Here ag^ is the RMS of counts of galaxies 




Figure 10. The halo model prediction for the large-scale galaxy 
bias bg L as a function of the averaged halo mass (Mjiost) 
^ = 0.1. Curves are for a different value of a (a = 0, 0.25, 0.5, 0.75 
and 1 from upper to lower). Open circles indicate a correspond- 
ing value for Mmin in intervals of A log Mmin = 1 staring from 
Mmin = 10^h~^MQ with Mmin = 10^^ Mq being marked by 
the filled circles. 

in spheres of 8/i^^Mpc radius and is derived from the best- 
fitting power-law model of the real-space correlation func- 
tion C = {ro/Sy via the relation {af^f = J2(7)(ro/8)^ with 
J2(7) = 72/[(3 - 7) (4 - 7) (6 - 7)2T] (Peebles 1980). Since 
af^ of the dark matter is a^^ ~ 0.9 at z — 0.1 (the mean 
redshift of the galaxy samples considered here) for as — 0.9 
(note that as is the linearly extrapolated value at z = 0), the 
corresponding bias is found to be bg,L — 1.2±0.1. Following 
the same procedure, we estimate a^^ of galaxy samples from 
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the SDSS. Adopting the best-fitting power-law models given 
in Table 2 of Zehavi et al. (2002), we find a^^ = 1.2±0.03 for 
their both "red" (defined by m — r > 1.8) and "high concen- 
tration" (c — rtjo/rso > 2.7) samples, giving hg^L — 1.3±0.03 
for z = 0.1 (the mean redshift of the samples). The lumi- 
nosity range of those galaxies is —22 < Mr — 5 log ^ < — 19 
with Mr,* — —20.8. Comparing those bias values (i.e., 
bg,L = 1.1 — 1.33) with the halo model predictions plot- 
ted in FigurJTUl we find (Mhaio) = (2 - 6) x IO^/i^^A/q for 
0-25 < a < 0.75 (Note that as is shown FigureCni the re- 
lation between the large-scale bias and (Mhaio) only weakly 
depends on a). This is broadly consistent with the EPS pre- 
dictions for the descendants of the OPEGs (with z' < 23.8 ) 
and also with that of the LBGs (with i' < 26.0). Therefore, 
we may conclude that, in the viewpoint of the mass evolu- 
tion of hosting haloes in the CDM model, the OPEGs and 
the bright LBGs are consistent with being the progenitor of 
the present-day early-type galaxies. 



5 SUMMARY AND DISCUSSIONS 

We have analyzed the two photometrically-selected galaxy 
samples from the deep mult-band images of the SXDS: the 
LBGs at z ~ 4 and OPEGs at 2: ~ 1. The contiguous large 
survey area of the SXDS enables us to measure the an- 
gular two-point correlation functions over a wide range of 
separations with a good statistical quality. Comparing the 
large-scale clustering amplitude with the corresponding halo 
model predictions, we have estimated a characteristic mass 
of hosting haloes. Then, adopting the EPS model, we have 
computed the predictions for the mass evolution of the host- 
ing haloes to explore the likely descendants of those galaxy 
samples. In particular, we have compared the predicted halo 
masses of these two populations at different epochs to make 
a possible evolutionary link. We have also examined ex- 
pected halo masses of the present-day descendants of those 
two galaxy populations. 

Our major findings are summarized as follows: 

(i) The measured bias functions (defined by equation |SJ| 
of both the OPEGs and LBGs consist of two parts: The 
decreasing small-scale part and the flat large-scale part. The 
transition scales between them are 1' and 10" for the OPEGs 
and LBGs, respectively. Those scales agree well with the 
virial radii of haloes with the characteristic mass of hosting 
haloes predicted by the halo model analysis in The shape 
of the bias function is thus understood by the standard halo 
model picture: The small-scale clustering arises from galaxy 
pairs located in the same halo, and the large-scale clustering 
arises from galaxy pairs located in two different haloes. 

(ii) The conclusive measurement of the fiat large-scale 
bias allows us to safely define the large-scale bias factor. 
Also, the large number of galaxies allows us to examine the 
luminosity dependence of the large-scale bias factor. It is 
found that the large-scale bias factors of both OPEGs and 
LBGs positively correlates with the luminosity (see bottom 
panels of Figures |21 and |3J|. 

(iii) Comparing the measured large-scale bias factors of 
the two galaxy samples with the halo model predictions, 
we estimate the characteristic mass of hosting haloes. The 
merit of using the large-scale bias factor is two- fold: First 
is concerning the measurement; since many bins are used to 



compute it, the measurement is less sensitive to a statistical 
noise in each bin. Second is concerning the halo model; the 
large-scale bias does not depend on Afi , and in addition, the 
resultant (Mhaio) is less sensitive to a. Therefore it provides 
with a reliable estimate of (A/haio) . The predicted character- 
istic halo masses of both galaxy samples are found to be pos- 
itively correlated with the luminosity. For the OPEGs, this 
may arise from a correlation between the stellar mass and 
the hosting halo mass, because the observed z'-hand mag- 
nitude (rest-frame _B-band) is well correlated to the stellar 
mass for population like OPEGs with little star formation 
activity. On the other hand, for the LBGs, it may suggest 
a correlation between the hosting halo mass and the star 
formation activity, rather than the stellar mass. This is be- 
cause the observed i'-band corresponds to ~1500A in the 
rest-frame, where the luminosity is most sensitive to the 
star formation, and in addition, the LBGs are generally in 
an active star formation phase. 

(iv) Utilizing the EPS model, we compute the predictions 
for the halo mass distribution of the LBGs' descendants at 
2 = 1 in the CDM cosmology, and then we compare it with 
the halo model prediction for the characteristic halo mass of 
the OPEGs. It is found that, for our fiducial subsamples (the 
OPEGs with z' < 23.8 and the LBGs with i' < 27.0), the 
typical hosting halo mass of LBGs' descendants is slightly 
smaller than the predicted mass of the OPEGs' hosting 
haloes. It is also found that the brighter LBG subsample 
(with i' < 26.0) is likely to evolve into the systems with 
halo mass compatible to the predicted one of the OPEGs. 
Therefore, we may conclude that, in the viewpoint of the 
mass evolution of hosting haloes in the framework of the 
CDM model, the bright (i' < 4- 1) LBGs are consistent 
with being the progenitor of the OPEGs. Accordingly, it 
seems less likely that the LBGs population, as a whole, has 
evolved into the OPEC population. 

(v) We also compute predictions for halo masses of the 
present-day descendants of both the galaxy samples using 
the EPS model. It is found that the predicted mass range 
for the LBG sample with i' < 27.0 is slightly but system- 
atically smaller than that of the OPEGs (with 2' < 23.8). 
On the other hand, the prediction for the LBG sample with 
i' < 26.0 agrees better with that of the OPEG descendants. 
In the latter case, the peaks of the PDFs are located at 
~ 2 X lO^^h'^MQ and the tail of the PDFs extends to the 
mass range of Af > 10^"^ h~ ^ Mq. Thus the present-day de- 
scendants of the bright LBGs and the OPEGs are likely 
to be located in massive systems such like groups of galax- 
ies or clusters of galaxies. We also estimate the character- 
istic halo mass of local early-type galaxy samples from the 
2dF and SDSS with the halo model, and it turns out that 
the predicted mass is in good agreement with the EPS pre- 
dictions for the present-day descendant's mass of both the 
bright LBGs and OPEGs. Therefore, it is concluded that, 
in the viewpoint of the mass evolution of hosting haloes in 
the CDM model, the OPEGs and bright LBGs are consis- 
tent with being the progenitor of the present-day early-type 
galaxies. 

One of the most interesting implications from the above 
findings is that the possible halo mass dependence of the 
LBG's star formation history. This is speculated from the 
above finding (iv) that the predicted descendant's halo mass 



© 2005 RAS, MNRAS OOO.ITHTni 



10 T. Hamana et al. 



of bright LBG subsample {i' < + 1) is found to be in 
very good agreement with the characteristic halo mass of 
the OPEGs, whereas it seems less likely that the faint LBG 
population has evolved into the OPEG population. The na- 
ture of the descendants of the faint LBGs at z ~ 1 is not 
clear but is naively expected that they evolve into other 
populations than old passively evolving galaxies with lower 
mass and bluer spectra. It is important to notice that a little 
star formation at z < 2 would be enough to push the galaxy 
colour outside of the OPEG's colour criteria (Yamada et al. 
2005). Therefore, the halo mass dependence of the epoch 
of truncation in star formation activity is one possibility 
of interpreting the finding (iv). It is also found that the 
UV luminosity of the LBGs correlates with the hosting halo 
mass [the above finding (iii)]. This may be an additional evi- 
dence for the mass dependent star formation history, though 
a connection between the above two findings, (iii) and (iv), 
is not clear. A possible scenario is that LBGs in more mas- 
sive haloes have more active star formation at significantly 
high redshift such as z ~ 4, and they turn into the pas- 
sive evolution phase earlier. This scenario can be tested by 
performing the same clustering analysis as presented in this 
paper progressively toward lower (and higher) redshifts. We 
might be able to see the transition at some point where star 
formation activities in massive haloes are truncated and the 
mass of the haloes that are hosting active star formation is 
being shifted to lower mass as time progresses. At the same 
time, on the theoretical side, we should explore a possible 
physical mechanism that drives such mass dependent star 
formation histories. 

Before closing, it is important to note that we have not 
argued that LBGs are the only path to become OPEGs 
at 2; = 1 or the present-day early-type galaxies. In fact, 
it is likely that the ancestor-descendant connection is not a 
one-to-one correspondence, and some different high-z galaxy 
populations may have evolved into a similar low- 2: popula- 
tion, and vice versa. For example, among the ancestors of 
the OPEGs at 2 = 1, there may be some haloes which have 
not coUpased by z = 4 and would be seen as LBGs at some 
later epochs, as well as the objects which are not UV lumi- 
nous enough to be selected as LBGs at z — 4 due to large 
amount of dust extinction and/or older stellar ages. This is 
one reason why we have not take into account the number 
density of our galaxy samples when we examine their pos- 
sible connection. We note, however, the number density of 
the bright LBGs nLBG(«' < 26.0) ~ 3.8 x IQ-^/i^'Mpc-^ is 
in fact comparable to that of OPEGs nopEciz' < 23.8) ~ 
4.7xlO"^/i^Mpc"^. Therefore, we could even argue from this 
comparison that the bright LBG-OPEG connection could 
be the major ancestor-descendant relation if the number 
density of bright LBGs would not decrease significantly by 
mergers. 
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